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We discuss photogenerated midgap states of a one-dimensional (ID) dimerized Mott insula- 
tor, potassium-tetracyanoquinodimethane (K-TCNQ). Two types of phonon modes are taken 
into account: intermolecular and intramolecular vibrations. We treat these phonon modes adi- 
abatically and analyze a theoretical model by using the density-matrix renormalization group 
(DMRG). Our numerical results demonstrate that the intermolecular lattice distortion is nec- 
essary to reproduce the photoinduced midgap absorption in K-TCNQ. We find two types of 
midgap states. One is a usual polaronic state characterized by a localized elementary excitation. 
The other is superposition of two types of excitations, a doped-carrier state and a triplet-dimer 
state, which can be generally observed in ID dimerized Mott insulators, not limited to K- 
TCNQ. 

KEYWORDS: strongly correlated electron systems, electron-phonon coupling, photoinduced phase transi- 
tion, polaron, density matrix renormalization group 



1. Introduction 

Midgap states observed in one-dimensional (ID) or- 
ganic band insulators are due to excited states that con- 
sist of localized carriers coupled with phonons.^'^ A po- 
laron state is a typical example of these midgap states 
and has attracted a lot of research interest. In particular, 
photogeneration of the polaron state has been intensively 
studied for several decades in the context of ultrafast gen- 
eration of a polaron and its relaxation dynamics. 

To understand the nature of these midgap states also 
provides us a good insight into the fundamental mech- 
anism of macroscopic photoinduced phenomena, called 
photoinduced phase transitions.'^'^ Currently, consider- 
able attention has been paid to the photoinduced phase 
transitions, in particular, of strongly correlated electron 
systems covering inorganic^^* and organic^^^'' materials 
including those which show photoinduced inverse spin- 
Peierls transitions. Although the transitions intrin- 
sically accompany large-scale changes, their initial seeds 
would be local and microscopic ones'^' ^ that are likely to 
be observed through the appearance of midgap states im- 
mediately after photoirradiation. Therefore, to examine 
the midgap states can help us to understand the starting 
point of the photoinduced phase transitions.^^ 

In this study wc focus on an organic material, 
potassium-tetracyanoquinodimethane (K-TCNQ). This 
material is classified as a ID dimerized Mott insula- 
tor below the spin-Peierls transition temperature, Tgp = 
395K.^^"^^ In 1991, Koshihara et al. have demonstrated 
that an irradiation of a pulsed laser weakens the lattice 
dimerization of K-TCNQ in the dimerized phase, which 
is regarded as a photoinduced inverse spin-Peierls tran- 
sition.^^ More recently, a high- resolution experiment has 
revealed several ultra-fast phenomena immediately after 



the photoirradiation, for example, the emergence of a 
midgap peak around 0.3 eV in the reflectivity spectrum 
within the time resolution of ISOfs."'^^' The decrease 
of the lattice dimerization is found to occur after the 
midgap state appears. 

The purpose of our study is to elucidate the nature 
of the midgap state in K-TCNQ. Our preceding study 
has made a suggestion that the midgap state is due to a 
purely electronic excitation of photoinduced mobile car- 
riers. This scenario is excellent in that the instant ap- 
pearance of the midgap peak is explained naturally. How- 
ever it is inconsistent with the experimental fact that 
the Drude component caused by the carriers does not 
appear. 

In this paper, we discuss an alternative mechanism; 
polaronic localized midgap states generated by the pho- 
toexcitation. We take account of two types of phonons, 
intermolecular and intramolecular vibrations, and treat 
these phonon modes adiabatically to analyze a theo- 
retical model by using the density-matrix renormaliza- 
tion group (DMRG).^^"^^ Our numerical results demon- 
strate that the midgap state in K-TCNQ can be re- 
produced by two polaronic lattice configurations. One 
is a configuration with the relaxed inter-molecular vi- 
bration mode and the unrelaxed intramolecular vibra- 
tion mode. This configuration provides a shoulder struc- 
ture around O.leV in addition to the main midgap peak 
at 0.3 eV. The other configuration is accompanied with 
the relaxed intermolecular and intramolecular vibration 
modes, where the shoulder and the peak merge into a 
single peak structure at 0.3 eV. We clarify the origins of 
these peak/shoulder structures and find out that spin ex- 
citations play an important role, which is characteristic 
to ID dimerized Mott insulators. 
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2. Model and Method 

2.1 Model 

In this work, we use the ID extended Hubbard model 
with Peierls and Holstein types of electron-phonon (e-ph) 
couphngs. The Hamiltonian is given by 

1=0 1=0 
N-2 N-1 

1=0 1=0 
J, N-2 J, N-1 N-2 

1=0 1=0 1=0 

where 

= J2^^Ul.<TCl,a + 4^Q + i,^), (2) 

and cj ^ (q,o-) creates (annihilates) an electron with spin 
a in the LUMO of a TCNQ molecule on site I, and 
ni.a = c\^ci_a. The Coulomb repulsion is taken into ac- 
count up to the nearest neighbor: the on-site Coulomb 
interaction U and the nearest-neighbor interaction V . We 
impose the open boundary condition (OBC) and N de- 
notes the system size. The parameter to gives the transfer 
integral of the uniform lattice. The distortion of the Z-th 
bond is denoted by ui, and the molecular deformation is 
given by vi. The e-ph coupling constants arc given by a 
and /3, and the elastic spring constants are denoted by 
Ka and Kp, respectively. Taking the adiabatic approxi- 
mation for phonons, we omit their kinetic energy terms, 
r is introduced to keep the chain length the same as that 
of the unrelaxed configuration, i.e., 

N-2 

E = 0- (3) 

(=0 

For simplicity, we re-define phonon variables in the fol- 
lowing; aui ui and (3vi —>■ vi. Accordingly, the elas- 
tic constants are renormalized as Ka/ct^ — > Ka and 

2.2 Method 

To obtain the stable (i.e., fully relaxed) lattice config- 
uration, we use the Hellmann-Feynman theorem; 

= and/or = for aU I. (4) 

dui ovi 

These equations lead to the following relations, 

ui = ((f,,,+i) +r)/X„, and (5) 

V, = ((n,) - l)/^/3, (6) 

where () is the expectation value for the eigenstate we 
obtained, and F is derived as 

N-2 

r = -E(*M+i>/(A^-i) (7) 

1=0 

from eq. (3). We evaluate the RHSs of eqs. (5), (6), and 
(7) by using the DMRG method for a certain lattice con- 
figuration (m, v) and F, and then we obtain new values 
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of (m, v) and F. The new configuration and F are used 
to evaluate the RHSs of eqs. (5), (6), and (7) again and 
we obtain next values. We execute this iterative proce- 
dure and finally reach the stable lattice configuration. For 
the final configuration, we apply the dynamical DMRG 
(DDMRG)^^"^^ to calculate the spectral function given 
by30 

^('^) = ^ E \mJ\^o)?K^ + Eo- E,), (8) 

i 

where 

J = i ^(to - - cl^Cl+l^a) (9) 

Lcr 

is the current operator, \ipi) the i-th eigenstate, and Ei 
the corresponding energy, xi^) is related to the optical 
conductivity spectrum as xi'^) = wcr(u;). The lattice con- 
stant is set to be unity for simplicity. In the actual com- 
putation, Lorentzians with the broadening parameter e 
are substituted for the 5 functions in eq. (8). We have 
carried out almost all DMRG calculations on = 40 
systems with retained bases m — 80. 

3. Optical Conductivity at Half Filling 

Before discussing polaronic midgap states in doped 
systems, we carry out calculation at half filling to obtain 
model parameters, U,V, and Ka for K-TCNQ. Kp is not 
determined in this study, as noted in §4. Our method of 
obtaining these parameters is as follows. First, we find 
the most plausible values of U and V by comparing the 
calculated optical conductivity ct^lu) with the experimen- 
tal result by Yakushi et al.'^^ Here the bond distortion is 
fixed to ui — {—lySto, where S is the dimerization pa- 
rameter and estimated a.s S = 0.46(= Se) from the ex- 
tended Hiickel calculation by K. Ikegami et al.^^ Thus, 
at this stage, we do not use the iterative procedure de- 
scribed in §2.2. 




Fig. 1. Optical conductivity spectra of K-TCNQ. The symbols 
show experimental result at 27 K in rcf. 31, and the solid line is 
our DDMRG result for TV = 40. 

Our best DDMRG result of (j{uj) and the experimen- 
tal data by Yakushi et al.^^ are shown in Fig. 1. The 
DDMRG resuh is calculated for U/to = 5, V/to = 1.7, 
and to — 0.291eV, with the broadening parameter e = 
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O.lBtQ. As for the on-site Coulomb interaction U, our es- 
timation is comparable to the values in the previous stud- 
ies. ^^"■^^ In addition, the relatively large V is consistent 
with our previous conclusion that large V is necessary to 
reproduce the sharp lowest peak and the broad higher- 
energy shoulder, '^^ which are characteristic to K-TCNQ. 

Then we determine Ka realizing S ~ 6c &s the stable 
lattice configuration at half filling to be Ka — 1.6/to- 
Actually, the obtained stable lattice configuration at half 
filling with OBC is not completely uniformly dimerized. 
Besides a small boundary effect due to the OBC, there 
is a small difference between even and odd bonds: the 
number of even bonds is larger than that of odd bonds 
by one. Hence, we obtain the relaxed lattice configuration 
where the average of dimerization is equal to dc] 



1 



N-1 



N-2 

E 

/=0 



(10) 



and then we confirm that (j{lo) for this stable configura- 
tion is almost identical to that for the uniformly dimer- 
ized configuration. 

4. Polaronic States 

Now, we discuss photogenerated polaronic states. In 
general, photoexcitation with energy corresponding to 
the charge transfer (CT) excitation generates two carri- 
ers. For ID Mott insulators, these are known as a holon 
and a doublon."^^ In the experiment by H. Okamoto et 
ai, the pump excitation with 1.55eV is considered to 
generate mobile holons and doublons.^'' These carriers 
move freely and are made separate within the time scale 
of h/to ~ 10^^"' sec, which is much shorter than the ex- 
perimental time resolution (~ 150fs)^^ and the phonon 
dynamics. Hence, we consider that the lattice relaxation 
starts after the carriers are well-separated, suggesting the 
interaction effect between the carriers is negligible. Be- 
cause of the particle-hole symmetry of this model, we 
consider the case where only one doublon exists on the 
system. We confirm that the case with a single holon 
shows the same results. As for the type of lattice relax- 
ation, we consider the following cases: {case 1) u-relaxed 
and w- fixed, {case 2) u- and w-relaxed, and {case 3) u- 
fixed and u-relaxed, where M(w)-relaxed means that the 
intermolecular (intramolecular) distortion u{v) is relaxed 
to the stable inhomogeneous configuration, and the u{v)- 
fixed denotes that the lattice configuration is fixed to the 
homogeneous one at half filling. 

Figure 2 shows uja{uj) of the polaron states for these 
three lattice configurations. The lowest-energy peak due 
to the open boundary condition is invisible in this and 
following figures. The results suggest that both of case 
1 and case 2 well reproduce the midgap peak at 0.3 eV, 
although a small difference exists between them. The for- 
mer shows a shoulder structure on the lower energy side 
of the main midgap peak. By contrast, case 2 seems to 
exhibit a single peak at 0.3eV. In case 3, there appears 
a midgap state above 0.5eV, being inconsistent with the 
experimental result. These observations suggest that the 
relaxation of u is necessary to reproduce the experimen- 
tal result. 
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Fig. 2. (Color online) Spectral functions for polaronic states 
where a single electron is introduced to Af = 40 systems. In 
cases 2 and 3, Kfs is set to 1.0/tQ. The inset shows the magnified 
picture for the midgap region. 



The remaining problem is whether the relaxation of v 
is necessary or not. In fact, we cannot reach the clear 
answer at present. The experimental resolution is not 
high enough to decide whether or not the low-energy 
shoulder is realized and to determine the plausible value 
of Kfj. In the following, we only examine the origins of the 
characteristic peak structures of case 1 and case 2, which 
are the candidates for the polaron state of K-TCNQ. 

5. Origins of Midgap States 

5.1 Case 1 

In Fig. 3, we plot oja{uj) in case 1, and that of the 
unrelaxed lattice configuration {case 0), where u and v 
are the same as those at half filling. Further we show 
U!a{Lu) for several intermediate lattice configurations be- 
tween case and case 1, defined as 

u{x) — u(case 0) + a;[M(case 1) — M(case 0)], (11) 

where the degree of relaxation is denoted by x, and in 
Fig. 3, X is shown in percentage. 
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Fig. 3. (Color online) Spectral functions for several lattice con- 
figurations (see the text). 



The most intriguing is that a midgap peak, termed 
(A), appears at 0.5eV in case 0. Since no analog is present 
in noninteracting band insulators, this midgap peak is 
characteristic to the ID dimerized Mott insulators. The 



4 J. Phys. Soc. Jpn. 



Full Paper 



Nobuya Maeshima 



peak (A) shifts toward lower energy as the lattice config- 
uration approaches the completely w-relaxed one (100%). 
We also note that the other midgap peak, termed (B), ap- 
pears in the lowest energy region for the slightly relaxed 
configuration (x = 20%). The peak (B) moves towards 
higher energy as the relaxation proceeds, and finally be- 
comes a shoulder at m(100%). 

To clarify the origins of these midgap peaks, we cal- 
culate the optical conductivity spectra of unrelaxed sys- 
tems with varying U and V, which are shown in Fig. 4. 
We can see three different classes of midgap peaks. One 
is the metallic low-energy component in the lowest fre- 
quency region. This is the Drude "precursor" in the OBC 
case.^^ Another excitation is observed around uj ~ 3to, 
which is well separated from the CT band for U/to > 8. 
This peak is caused by the intradimer CT excitation at 
the carrier site, while the carrier is mobile. In ref. 25, 
we argued that this CT excitation is a possible origin of 
the midgap peak in K-TCNQ. However, the present es- 
timation of the model parameters for K-TCNQ {U/tQ = 
5, V/^o — 1-7) demonstrates that the corresponding peak 
is hidden by the CT excitations around uj ^ AtQ. 
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Fig. 4. (Color online) Optical conductivity spectra for several sets 
of (U,V) with unrelaxed lattice configuration (case 0). The re- 
sults are obtained with Lanczos diagonalization for A'^ = 12 sys- 
tems with 7 electrons and the OBC imposed. The broadening 
parameter is set to e = O.lto. The filled circles show the singlet- 
triplet gap given by eq. (17) 

Now we focus on the other class around oj/to = 1 — 2. 
For the parameter set corresponding to K-TCNQ, this 
peak is located at uj/to ^ 2, which is identical to the 
midgap peak (A) . Thus this midgap peak (A) can appear 
generally in ID dimerized Mott insulators. It should be 
noted that its location is quite sensitive to U and V; 
the excitation energy increases as U decreases and as V 
increases, implying that the midgap peak (A) is caused 
by spin excitations. 

In fact, it can be shown that the excitation caus- 
ing the peak (A) includes a singlet-triplet excitation. 
To show this, we use the analytical expression in the 
decoupled-dimer limit. In the 0-th order approxima- 
tion, the ground state of the system is represented by 
the direct product of the singlet-dimer states; 

|7Ag) = |G)o®|G)i®---®|G)^,_i, (12) 



where \G)i, {I — 0, - ■ ■ , Nd ~ 1) is the singlet-dimer state 
on the ^-th dimer, and Nd = N/2. Then the one-electron- 
doped system is denoted by the linear combination of a 
direct product given by 

\G)o |G)i • • • |0^)„ (g) • • • |G)Ar,_i. (13) 

Here, |0|)„ has one more electron, and is defined as 

(14) 

with spin a. Then we can construct a CT excitation via 
a interdimer CT process to obtain 

|G)o |T)„ ® |Oj)„+i • • • |G)jv,_i, (15) 

where 

|T^)n = 4T4+iTlO)- (16) 
A schematic picture of this excitation is given in Fig. 5 
(a). Comparing eq. (13) and eq. (15), we find that the 
corresponding excitation energy A is equal to the singlet- 
triplet gap given by 

A = + ^{U^Vy/A + itlil + Sy. (17) 

In Fig. 4, we plot A for each parameter set and find that 
A explains the excitation energy of the midgap peak (A). 




Fig. 5. Schematic energy diagram for relevant midgap states (a) 
and (b) causing peaks (A) and (B), respectively. 

The existence of two states, |r)„ and |0|)„i, is con- 
firmed in the following analysis. If the state causing the 
peak (A) has these states, the doubly occupied site and 
the up spin must be separated (n ^ m). This is com- 
pletely different from the states that include only |0|)n, 
where the doubly occupied site and the up spin are lo- 
cated on the same dimer n. To demonstrate this spin- 
charge "separation" , we calculate the corresponding cor- 
rection vector, ■^^ 

^ ^ ^ ■^iV'o), (18) 

with normalization constant A. Figure 6 shows the 
charge and spin density profile of the correction vec- 
tors for a nearly completely relaxed configuration (80%). 
Here, we used the broadening parameter e — 0.16fo, 
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whose result is found to be almost identical to that with 
e = O.OSto (not shown). We can see that the spin density 
of the state |0.32eV) concentrates on the center of the 
lattice deformation and the charge density has a double 
peak beside the center, which shows that the spin and 
charge densities are separated as discussed above. 




Fig. 6. (Color online) Charge and spin distributions in several 
states for x = 80% with N = 38. 



By contrast, the spin and charge densities of the state 
|0.12eV) have similar distributions; both have double 
peaks and nodes at the center. Since the lowest state 
|0) has no node, the state causing the midgap peak (B) 
is the second lowest state made of the elementally exci- 
tation |0|), as shown in Fig. 5 (b). 

5.2 Case 2 

Now, we turn our attention to case 2. Figure 7 shows 
the spectral function Lua{uj) in case 2 for several values 
of K/3. We can see that the shoulder structure (B) for 
Kpto = oo approaches the main peak (A) as Kfj de- 
creases. Finally, these seem to merge into a single peak. 
This is because the excitation gap between the lowest po- 
laron state and the second lowest state increases as the 
one-electron potential due to the intramolecular defor- 
mation V becomes deep, as shown in Fig. 8. By contrast, 
the intermolecular deformation u changes slightly, keep- 
ing the singlet-triplet gap A almost constant. 
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Fig. 7. (Color online) Spectral functions in case 2. 




Fig. 8. (Color online) Lattice configurations for case 1 and case 
2 with Kpto = 1.0 for Af = 38 systems. 



6. Summary and Discussion 

We have studied polaronic midgap states in the ID 
extended Hubbard model with Peierls and Holstein types 
of electron-phonon couplings to clarify the origin of the 
midgap peak in K-TCNQ. Our main conclusion is that 
the relaxation of the intermolecular lattice distortion u 
is necessary to reproduce the observed midgap peak at 
0.3 eV. 

The midgap peak is numerically reproduced in the 
spectral function uja{w) for two polaronic lattice configu- 
rations. One has relaxed u and unrelaxed intramolecular 
vibration mode v {case 1), and the other has relaxed 
u and relaxed v {case 2). The former configuration is 
found to exhibit an additional shoulder structure on the 
lower energy side of the midgap peak. In this case, we 
have clarified the origins of the shoulder and the midgap 
peak. 

The shoulder structure is caused by the excitation from 
the lowest bound state to the second-lowest state, both of 
which consist of the elementary excitation |0|). This is 
similar to the polaron-based lower midgap peak of usual 
ID band insulators, where the peak stems from the ex- 
citation from the lowest bound state to a higher state of 
introduced carriers (holes or electrons).^' 

By contrast, the midgap peak is characteristic to ID 
dimerized Mott insulators. The peak corresponds to the 
interdimer CT excitation from the state |0|) to a super- 
position of the state |Oj) and the triplet dimer state |T). 
Its excitation energy is shown to exhibit spin-excitation- 
like (t/, y)-dependence. 

Our results provide some suggestions to the photoin- 
duced inverse spin-Peierls "transition." First, the 
polaronic lattice relaxation proceeds within the exper- 
imental time resolution (=150fs). This implies that the 
phonon motion to generate the polaronic state is much 
faster than the observed coherent oscilation.'^^ To dis- 
cuss the speed of the phonon motion, we here evaluate 
the bare phonon frequency. The experimentally observed 
dimerization length (~ 0.165A )'^^ and the dimerization 
strength of the transfer integral (^^o — 0-46 x 0.291eV) 
lead to the e-ph coupHng constant a ^ 0.81 cV/A. 
Then the unrenormalized elastic constant Ka is obtained 
to be 3.6 eV/A^. These values for physical parameters 
give the bare phonon frequency lo — 2 Ka / Mtcnq ^ 
2.6 X lO^^s-^ where Mtcnq ~ 3.4 x lO^^^kg is the mass 
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of a TCNQ molecule. Thus the speed of the bare phonon 
sub ps) is found to be comparable to the experimental 
resolution time. It remains unresolved whether the bare 
phonon frequency determines the speed of the polaron 
formation or not. 

Second, the photo-induced melting of the dimeriza- 
tion is explained by the "impurity (polaron) doping" pic- 
ture, not by the mobile-carrier picture, where the mo- 
bile photodoped carriers destroy the singlet dimerization 
of the (iimcrized phase. Third, one carrier can weaken the 
lattice dimerization roughly over 10 molecules as shown 
in Fig. 8. Assuming that interaction between doped carri- 
ers is negligible, we estimate that the dimerization order 
decreases linearly with respect to doping concentration 
at least up to 0.1 carrier per molecule. This is equivalent 
to 0.05 photon per molecule, which is equal to the pho- 
ton density where the AR/R at 0.71 eV saturates [see 
Fig. 2(d) of ref. 17]. 

Our estimation of the model parameters for K-TCNQ 
has demonstrated that K-TCNQ does not show the ab- 
sorption peak caused by the conversion of photogen- 
erated elementary excitations, which was proposed as 
the origin of the midgap peak.^^ This absorption for K- 
TCNQ is hidden by the CT excitation because the CT 
gap is nearly equal to the location of this absorption. 

The existence of this conversion can be confirmed in 
other quasi-lD materials that have a larger CT gap. For 
example, an organic radical crystal, l,3,5-trithia-2,4,6- 
triazapentalenyl (TTTA) is a candidate of these materi- 
als. In ref. 34, we have estimated model parameters for 
TTTA as 5 = 0.4 and U/to = 10, V/to = 1.7, to = 0.235 
eV. Since this material has relatively large U/to, the ab- 
sorption caused by the conversion can be well separated 
from the CT band. 
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